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FORWARD 


This  report  presents  a  comparison  between  a  simplified  equation 
for  stagnation  point  heat  transfer  and  the  experimental  results  of  an 
investigation  on  heating  rates  to  cylinders  and  spheres  subjected  to 
flame  impingement  from  a  l/10-scale  model  of  an  Atlas  vernier  engine, 
carried  out  under  contract  No.  AF-OU(6li5)-h. 
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NOME  FUTURE 

A  s  cross  sectional  are a  of  calorimeter 

c  -  specific  heat  of  material 

*  •  » 

_  specific  heat  of  gases  at  constant  pressure 
D  -  body  di arc* ter 
De  diameter  of  nozzle  exit 
g  =  gravitational  constant 
li  =  enthalpy 

k  •  conductivity  of  material 

"fc  •  ratio  of  densities  across  a  normal  shock 

P  -  pressure 

q  -  heat  transfer  rate 

7?  ■  universal  gas  constant 

r  ■  radius  of  body 

rs  -  outside  radius  of  shell 

T(n  -  inside  radius  of  shell 

T  ■  temperature 

t-  -  thickness 

•  equivalent  thickness  defined  by  E^*.(37)or  (38) 

U e  "  velocity  at  the  outer  edge  of  the  boundary  layer 
w  -  weight  of  calorircoter 

x  ■  distance  along  centerline  measured  from  nozzle  exit  (Section  III) 

x  •  distance  along  surface  of  body  (Section  I) 
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Oreek  Symbols 

*  m  tyc  thermal  diffusivity  of  material 
ot(TS/Tj-  function  given  by  Eq.  (8) 
c/'  —  function  given  by  Eq„  (Hi) 

c*  "  -  function  given  by  Eq„  (17) 
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r  - 

G  - 
K  - 

/■  * 

t  - 

e  - 

*f  • 

oc  = 


£  . 
s  - 

t  - 

w  - 

oc  • 

* 


ratio  of  speci fic  heats 

angle  between  radius  vector  and  centerline  of  body 

constant  in  the  temperature-viscosity  relation 

viscosity  in  lb/sec«  ft 

angle  about  the  axis  of  body 

density  of  gas  or  material 

tine  in  sec. 


fraction  of  heat  losses  by  conduction  to  total  heat  input 
at  T\;  -  620°R 


ubscripts 


chamber 

conduction  losses 

at  the  stagnation  point,  also  stagnation  conditions,  unless  otherwise 
noted. 

total 

wall  conditions 

conditions  ahead  of  shock 

evaluated  at  a  reference  temperature 
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A  review  of  the  available  theoretical  expressions  on  the  velocity  gradient 
at  the  stagnation  point  of  a  cylinder  or  sphere,  (an  important  parameter 
in  stagnation  point  heat  transfer)  discloses  no  ap  reciable  discreoancy 
between  the  various  theories.  The  equation  of  Fay  and  Riddell,  derived  for 
the  case  of  dissociated  air,  is  next  transformed  by  using  the  perfect  gas 
law  and  the  available  data  on  the  viscosity-temperature  relation  of  the 
exhaust  gases  of  the  vernier  engino.  A  simplified  equation  giving  the  heat 
transfer  rate  as  a  function  of  the  body  diameter,  the  Mach  number  ahead 
of  the  shock,  and  the  stagnation  or  chamber  pressure  is  thus  obtained. 


As  a  test  of  the  validity  of  this  procedure,  the  incompressible  flow  equation 
given  by  Sibulkin  is  modified  by  the  incorporation  of  the  reference  tempera¬ 
ture  method  of  Eckert.  This  is  compared  with  the  previous  result*  The 
maximum  deviation  of  this  result  from  that  previously  determined  is  found 
to  be  less  than  U  percent. 

In  Section  II,  two  principal  sources  of  error  in  the  experimental  data, 
the  error  due  to  the  curvature  and  thickness  of  the  calorimeters  and  that 
due  to  conduction  losses,  are  discussed,  and  a  method  for  the  elimination  of 
these  errors  is  given.  It  is  also  shown  here  that  the  discrepancy  between 
theory  and  experiment,  with  regard  to  the  variation  of  the  heat  transfer 
rate  with  wall  temperature,  can  be  explained  easily  when  the  conduction  losses 
are  taken  into  account. 

In  Section  III  a  comparison  is  made  between  the  theoretical  and  experimental 
results.  It  is  found  that  the  adjusted  experimental  values  are  in  good 
agreement  with  the  theoretical  curve,  while  the  unadjusted  data  deviate 
appreciably  in  the  region  of  low  heat  transfer  (high  Mach  number)  of  the 
jet,  where  the  conduction  losses  become  important. 
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I  THEORETICAL  HEAT  TRANSFER  RATES 


1  -  Velocity  Gradient  at  the  Stagnation  Point  of  a  Blunt  Body 

The  velocity  gradient  at  the  outer  edge  of  the  boundary  layer  near 
the  stagnation  point  is  required  in  order  to  compute  the  heat  transfer 
rate  as  shown  in  Section  2  below0  A  number  of  available  theoretical 
treatments  give  approximately  the  same  value  for  the  velocity  PTadient0 

L.  Leos  and  Fay  &  Riddell  (Reference  1)  use  the  hypersonic  or  New¬ 
tonian  flow  approximation  which  leads  to  the  expression 


(&>. 


*  2-7 IF 

3>  v  P* 


d *  s 

Using  the  perfect  gas  law  this  equation  takes  the  form 


(1) 


(—1  = 


2-  \Jzax 


(la) 


Hayes  (as  reported  in  Reference  2)  gives  an  expression  for  the 
velocity  gradient  in  terms  of  1c  ,  the  density  ratio  across  the  normal  shock j 


(4^)  = 

'd%  s  D 


(2) 


Using  the  definition  of  the  speed  of  sound  and  expressing  £  in  terns  of  the 
free  stream  Mach  number  one  obtains,  after  some  manipulation. 


(dii)  =  z_, TW  \fz rts 


‘s  ~  v  r*'  ’  "  (2a) 

Finally,  Li  and  Oeiger  (Reference  2)  give  a  modified  expression 


(t%  TT  Sj *(»■-*> 


d%' s  D  v  '  05 

After  some  manipulation  this  equation  can  be  put  in  the  following  form: 


(r-*)  =  f  x/1*  r'- 

'  d+  's  I>  \  X+i  L 


t-i  _  \ 

2lf+i)  (tT+i) 


.]  V2RTS(3.) 
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Eqa  (2a)  and  (3a)  are  of  the  came  form  as  Eq.  (la)  except  for  a 
small  factor  which  forTf"  1.267  is  equal  to  1.055  in  the  case  of  Eq.  (2a), 
and  1«015  in  the  case  of  Eq.  (3a)  with  M-5.  These  equations  are  plotted 
in  Figure  1.  Since  the  velocity  gradient  appears  as  a  square  root  terra 
in  the  heat  transfer  equation,  the  discrepancy  in  using  any  one  of  the  above 
equations  will  be  small. 


2  .  Simplified  Heat  Transfer  Equation  for  a  Perfect  Qaa 

Fay  and  Riddell  (Ref.  1)  give  an  equation  for  heat  flux  to  the  stagna¬ 
tion  point  of  a  body  of  revolution  in  dissociated  air  which,  for  Prandtl 
number  equal  to  0.71,  is 

j  •  0.91*  (p-yU-)  (p,H,r  (h-h*)  <ll) 

Recently,  Beckwith  pointed  out  tha*  the  above  equation  also  c  orrelates  the 
results  of  perfect-gas  solutions,  provided  that  Sutherland  formula  for 
viscosity  is  used. 


Equation  (U)  takes  the  following  form,  after  multiplication  and  division 
by  (psuxf and  substitution  of  Eq.  (la) 

'  so  [e>“> 

Using  the  perfect,  gas  law  and  considering  that  the  pressure,  p,  is  constant 
across ,and  at  the  outer  edge  of  the  boundary  layer f one  obtains 


o.i 


^  (c^ 11  ~  *+* r~)(***Q (6) 

or,  grouping  all  the  terras  that  depend  on  the  teraperatures,  Ts  and  Tw  , 

(7) 


*(Tltr~) 

;  /D 


where 


*(T,T.)  =  1.33  ('£$  )  (Ch  Ts  -  C^Ti)  (2?  KZ) 


(8) 
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For  hypersonic  flow  the  value  of  pis  given  approximately  by 


whore  j from  the  one-dimensional  flow  relations 


A.  = 


h 


(  I  4  tl  H.‘ ) 


(9) 


(10) 


therefore 


t,  r 

•  /  J-  i  .  .  v  \tyf-l 


(l  +  tj!  Mi) 

Eq.  (7)  then  becomes,  after  substitution  of  Eq.  (11)  t 


(11) 


/9  iu. 


Q  _  ov  i/7a  ^ 
i  Wo  ( /+ 1!  3“'*“ 


(12) 


ortif  ^  is  given  in  lbs/in  ,  and  D,in  inches 


1-  - 


a 'STS 


where 


SB  (it  !M  Mi) 

•e  * 


0('  =./  492  >  <X 


/itu 

7>r  ftx 


(13) 


(Hi) 


For  )f*1.267  and  /?«.  «  3li7  psia,  q  becomes  a  function  of  Mach  number 
and  diameter  only. 


f  (  D  tmhfi J  (l  ) 


where 


'  **  (,+<,,»*  Mir1 

ci  •  3.13  x  106  «  ^ 

•  U.16  *  106  (r§.  fuTTL  (tfWn-CfJjtyW)  (16) 

In  a  slightly  different  form 

W'  "  4'T";  (17) 


-  t  t  «  -trt* 
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Eq.  (17)  it»  plotted  in  Fig.  2  for  a  range  of  value?  of  Tw  using  Tfl  as  a 

parameter,  and  the  available  values  of  viscosity  and  specific  heat  of  the 
exhrjst  gases.  It  is  noted  from  this  equation  that  the  effect  of  wall 
temporature  on  the  heat  transfer  rate  is  contained  in  the  temperature  gradient 
tern,  the  terra  )0#*  being  insensitive  to  a  variation  of  of  a 

few  hundred  degrees.  Thus,  in  the  range  of  T«  •  600  -900°R  and  Tg  •  5000  - 
6000°P#  Eq.  (17)  can  be  simplified  to  the  following  equation,  accurate  to 
+  1,$  percent! 

o(  -  3,65  x  10-*  ( )  (18) 

and  tho  heat  transfer  rate,  from  Eq.  (15),  becomes 

,  ^_«$47p»A* 

flip* 

r  l/i)  (UOJUS-M.)  />,„  !M*es  I 

For  Tg  •  5500°R,  Tw  •  620°P>  and  using  the  corresponding  values  for^u  and 
C  ,  one  finds  c<*-  9.55  x  10°  ,  so  that 


q  » 

/  VS 


f>±  -  3f  7  pi 

J  Ti  •  S5ac  *4? 
r„  «  czo 

J>  to  incAez 


The  above  equations  apply  to  the  case  of  a  sphere.  For  the  case  of  a 
cylinder,  the  equivalent  equations  are  obtained  by  dividing  by  )/2~, 


3  -  Comparison  with  the  Equation  Obtained  by  the  Reference  Temperature 
Method. 

It  has  been  suggested  that  the  incompressible  stagnation  point  heat 
transfer  equation  given  by  Sibulkin  can  be  used  for  compressible  flow  pro¬ 
vided  that  the  gas  properties  are  evaluated  at  a  reference  temperature 
given  by  Eckert.  (See  Ref.  3)»  Thus,  the  modified  equation  is 


U  -  ^1%  v'TP  (h.  -  D 


fr.r.  ^)OX  *  .  - 

which  f or  fr  ■  0.71  becomes 


fa.  =  ’0.9*S  l/f  V  (h,-h„)  V^/dx). 


(21a) 
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where  the  reference  temperature  ie  computed  from  the  simplified  equation 

T*  •  I  <  T«  ♦  I„  )  (22) 

Comparison  of  Eq*0  (21a)  and  (U)  gives 

frt  ^  \/  (23) 

1 

Using  the  perfect  gas  relation,  Eq.  (23)  may  be  modified 


as  follows: 

H  . 

!  •  ■  ? 

(t  DV# 

J  II 

(2U) 

y  ■«>  t- 

/  M***  <5*  )  °'1 

T«f 

A  temperature  -  viscosity  relation  of  the  form 


may  be  assumed,  and  Eq.  (2li)  becomes 


j 

(25) 


(26) 
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A  value  of  X-  0.62  to  0.68  was  obtained  from  the  viscosity  temperature 
data  of  the  exhaust  gases  in  the  range  of  Tw  ■  620  -  900°R*  Using  this 
value  of  Kin  Eq.  (26)  the  ratio  of  the  heat  transfer  rates  given  by  the 
two  equations  was  found  to  vary  slightly  from  unity 

9r*  /.o>s  t  /•  ose 

f 

Thus  Eqs.  (21a)  and  (U)  give  approximately  the  sane  heat  transfer  rate. 
Values  obtained  from  Eq0  (21)  have  been  plotted  in  Figure  (5). 
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II  ANALYSIS  OF  THE  EXPERIMENTAL  RESULTS 

1  -  Qeneral  Equation  of  Heat  Conduction  in  the  Calorimeter 

The  differential  equation  for  unsteatty-state  heat  conduction  in 
spherical  coordinates  is 

j7~  r  i  o>  (rx  dT)  .  j _  A  (&n  9  i^)  4. 1 —  ell Z~  1 

Jp  -  &  |  rZ  '  dr)+  r*  imp  4*  "  **  r*sm$  J  (27) 

Assuming  that  all  derivatives  of  T  with  respect  to  9  and  f  are  not  functions 
of  r,  integrating  with  respect  to  and  using  the  boundary  condition 

at  T"a  /'I*  one  obtains  the  temperature  gradient  at  the  outer  surface, 

r*rs  , 


/  £l)  ,  G  -  C,  £T  _  r*~r'"  ±  ( sin 9  £T) 

e>r  s  3<y/jl  si»&  06' 

The  heat  transfer  rate  is  then  given  by 

f/jl  _  _ r/-  't— "cT  .  ic'  J)_  /stn&J?T) 

/>L;  (7? { -  ~r  + 


The  corresponding  equation  for  heat  transfer  to  a  cylinder  is  found  in  a 
similar  manner  starting  from  the  differential  equation  in  cylindrical  coordi¬ 
nates 


±  A  (riT)  +  J-  £1] 
jr  '  l  r  dr'-  dr1  r<-  OB'  J 


Thus 


\al  .  -  -  pc  ll-TnT  £T  +  i- 

'  2  /;  ^  a  &  tr> 


The  first  ter*  on  the  right  hand  side  of  Eqs,  (29)  and  (31)  represents  the 
heat  storage  in  the, shell  element  and  the  second  term  »  the  conduction  away 
from  the  element.  When  the  site  of  the  element  is  taken  equal  to  the  size 
of  the  calorimeter  this  second  terra  represents  the  heat  losses  by  conduction. 
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2  -  Correction  for  the  Thickness  and  Curvature  of  the  Calorimeter 

■■  11  1  ■  —  —  —  r  m  m  . -  1  11  ’■ 1  1  ■  •  —  1  - -  - - 

The  equation  uaed  in  computing  the  heat  storage  in  the  calorimeter  from 
the  experimental  temperature-time  curve  isi 


r.  JT 

A  dr 


This  equation  applies  strictly  to  the  case  of  a  flat  plate#  For  a 
calorimeter  haring  the  form  of  a  spherical  shell  the  storage  term  in  Eq.  (29) 
applies  i 
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which  can  be  put  in  the  form 

q  =  pet  fcf  dS  v  c  AJ  (lx)  (35) 
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Comparison  of  Eq„  (35)  and  (32)  gives 

=  (jV  -a  (36) 

Eq,  (36)  applies  also  to  the  case  of  a  cylindrical  shell  element;  however 
for  a  sphere,  one  finds,  after  some  manipulation  of  Eq,  (3U), 

=  /  _  i?  *  J  f1-  (37) 

whereas,  for  a  cylinder,  using  the  storage  term  of  Eq.  (31)* 

U^J  1  =  /  -  Jl  (38) 
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To  compute  the  losses  by  conduction  a  knowledge  of  the  surface  tempera¬ 
ture  distribution  is  required.  Since  no  experimental  values  are  available, 
a  temperature  distribution  of  the  form 

7"«  CCi  &  (39) 

may  be  assumed.  This  distribution  checks  well  with  experimental  values 
obtained  by  the  NACA  at  a  Mach  number  of  six,  and  for  small  values  of 
&  (For  a  calorimeter  dioraeter  of  0o128B  and  a  sphere  diameter  of  3/b% 

&  is  approximately  9*6  degrees)  The  temperature  derivative,  appearing 
in  the  conduction  term  of  Eq.  (29)  is  found  from  Eq.  (39) 


-L  d.  (sm  &  dT )  ~  -  2  Xi  ccsfr  *  -  2 
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and  the  he?t  loss,  using  the  conduction  tern  from  Eq.  (29),  is 

/,  =  -Z  T„.J  (sfik.r*)  ()a) 

The  above  equation  and  the  temperature  distribution  given  by  Eq.  (39)  are 
strictly  valid  for  a  homogeneous  sphere.  In  the  case  of  a  calorimeter 
material  different  from  the  model  material  Eq.  (bl)  can  be  used  for  esti¬ 
mating  purposes.  With  the  value  of  k  for  the  sphere  material.  The  validity 
of  this  assumption  is  demonstrated  by  an  indirect  approach  in  Section  b 
below. 

For  tho  case  of  the  cylinder. 
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U  -  Effect  of  Conduction  Losses  on  the  Heat  Transfer  Rates  at  Different 
Wall  Temperatures! 

The  large  discrepancy  between  experimental  results  and  theory  with 
respect  to  the  variation  of  the  heat  transfer  rate  with  temperature,  shown 
in  Fig.  1,  page  5  of  Reference  b,  can  bo  accounted  for  by  including  the 
heat  losses  by  conduction.  The  total  heat  transfer  from  the  gas  to  the 
sphere  at  the  stagnation  point,  for  wall  temperatures  of  62D°R  and  TWs  is 
respectively 

jit  J*  =  +  (y, 

%,  =  P‘f  ^  *  VT-'  (US) 

From  Eq.  (bl)  it  is  seen  that  the  heat  losses  are  proportional  to  the 
temperature  and  therefore 
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Moreover  from  the  theoretical  heat  transfer  results 
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Combining  Eqs.  (UU)  and  (b5)  substituting  Eqs.  (b6)  and  (h7)  and  rearrang¬ 
ing  one  obtains 
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Eq.  (U8)  plotted  in  Fig.  U  forX"  0,  0.10,  0.15 

To  ietermine  the  value  of  \  for  the  conditions  of  Fig.  1  of  the  above  c 
mentioned  report,  the  heat  loss  was  computed  from  Eq.  (Ul)  to  ber,0.71  x  10^ 
Btu/hr  ft^  and  v&s  found  from  theory  at  M>e*  11  to  be  b.6  x  ICk  so  that 
X  *  0.15.  The  experimental  curve  is  also  plotted  in  Fig.  It  as  a  dashed  line. 
It  coincides  with  the  theoretical  curve  for  X  ■  0.15,  which  is  the  same  as 
the  computed  value.  This  agreement  indicates  that  the  use  of  Eq.  (Ul)  for 
computing  losses  ie  justified.  The  use  of  the  conductivity  of  the  sphere 
material  (steel)  may  also  be  justified  from  the  fact  that  the  conductivity 
of  the  calorimeter  material  (copper)  is  much  larger. 
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III  CORRELATION  OF  THEORY  AND  EXPERIMENT 


The  experimental  results  reported  in  Reference  (7)  must  be  corrected  for 
the  errors  discussed  in  Sections  2  and  3.  Using  fy+cio  to  denote  the  data 
as  reported,  the  adjusted  heat  transfer  rate  is  given  by. 

fa  =  (fV  +  !“■>“  m 

For  a  3A  in-diameter  sphere  with  shell  thickness  t  *0.10  in.  this  equation 
becomes 


fa  -  &7(,  +  7 l.ooo  (50a) 

The  experimental  points,  thus  adjusted,  are  plotted  in  Fig.  5.  The  data 
on  models  with  insulated  calorimeters  has  been  adjusted  only  for  the  error 
due  to  curvature.  The  simplified  equation,  Eq.  (20)  in  Section  I,  has  also 
been  plotted  in  this  figure  using  the  variation  of  Mwith  VdJ, See  Ref.  5) 
shown  in  the  lower  portion  of  the  curve.  The  agreement  between  theory  and 
experiment  is  seen  to  be  satisfactory.  The  experimental  points  have  also 
been  plotted  in  Fig.  6  without  any  correction  and  compared  to  the  theoretical 
results.  It  is  seen  that  heat  losses  cause  an  appreciable  error  at  the 
high-Mach-number  portion  of  the  curve  where  the  heat  input  is  low.  The  test 
results  on  a  3/U"  cylinder  are  plotted  in  Figure  7  along  with  a  plot  of  values 
of  o  obtained  from  Eq.  (20)  divided  by  kT.  The  experimental  points  are 
plotted  as  reported  in  Reference  (6)  without  any  correction  for  the  curvature 
and  thickness  of  calorimeter  and/or  for  conduction  losses.  From  Fig.  3  *r/ 1 
is  approximately  equal  to  0.87,  for  a  cylindrical  shell  element.  However,  the 
calorimeter  used  wa3  circular  and  this  correction  cannot  be  made  with  cer¬ 
tainty.  Moreover  a  uniform  constant  was  used  for  all  calorimeters  in  reducing 
the  experimental  data,  even  though  this  constant  was  expected  to  vary  slightly 
with  each  calorimeter.  The  heat  loss  computed  from  Eq.  (U3)  at  a  wall 
temperature  of  620°R  is  1*1,000  Btu/h*1  ft  .  This  correction  will  raise  the 
values  of  the  experimental  points  at  high  Mach  numbers,  and  give  a  better 
agreement  with  the  theoretical  curve.  It  is  emphasized  that  corrections 
for  curvature  and  heat  conduction  have  not  been  applied  to  the  data  presented 
in  Fig.  7. 

Finally,  it  should  be  pointed  out  that,  in  view  of  the  appreciable 
rarefaction  of  the  Jet  at  a  Mach  number  of  10  or  higher  (See  Fig.  A-2  p.  1*6 
of  Reference  8),  slip-flow  effects  will  tend  to  reduce  the  heat  transfer  rates 
to  cylinders  and  spheres,  and  the  theoretical  curve  is  not  strictly  applicable 
beyond  this  Mach  number.  From  Fig.  D-l*  (p.  71)  of  this  reference  it  is 
estimated  that  there  will  be  a  1J>  percent  reduction  in  the  average  heat 
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transfer  rate  to  b  3/U"  sphere  at  M  ■  10.5.  Thu*,  the  experimental  results 
in  this  rarefied  region  should  be  expected  to  be  somewhat  lower  than  the 
theoretical  curve,which  was  based  on  continuum- flow  theory. 
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